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Abstract. Newtonian dynamical systems which accept the normal shift on an arbitrary Riemannian manifold are considered. 
For them the determinating equations making the weak normality condition are derived. The expansion for the algebra of 
tensor fields is constructed. 
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1. Introduction. 



The concept of dynamical systems accepting the normal shift was introduced in [1] (see also [2] 1 ). It appears as 
the result of transferring the classical Bonnet transformation from geometry to the field of dynamical systems. In [1] 
and [2] the Euclidean situation was studied i.e. the dynamical systems in R" accepting the normal shift (we refer the 
£NJ ' reader to that papers for the detailed bibliography). 

Apart from recent results of [1] and [2] one can find another purely geometric generalization of Bonnet transformation 
from [3] and [4] which is realized as a normal shift along the geodesies on some Riemannian manifold. The natural 
question here is how do these two generalizations relate each other. This question was investigated in [5]. There the 
| special subclass of dynamical systems accepting the normal shift was separated for which the normal shift is equivalent 

to the geometrical generalization of the Bonnet transformation for some conformally- Euclidean metric in R™. Such 
systems are called metrizable. Comparing the explicit description of metrizable dynamical systems given in [5] with 
the examples of [1] and [2] one can conclude that non-metrizable systems do exist. Therefore the the concept of 
normal shift along the dynamical system is wider than the Bonnet transformation for conformally-Euclidean metrics. 
However it do not embrace the case of Bonnet transformation for arbitrary metric. In this paper below we consider 
the most general situation and study the dynamical systems accepting the normal shift on an arbitrary Riemannian 
manifold. 



2. Newtonian dynamical systems on the Riemannian manifold. 
The main object considered in [1] and [2] is the second order dynamical system in R™ 

(2.1) f = F(r,r) 

reflecting the Newton's second law. Now let r be not radius- vector of a point in 1" but the vector of local coordinates 
for some manifold M. This case let us write the equations (2.1) in form of the system of the first order equations 

(2.2) x i =v i i) l = $ l (r,v) 

Systems of differential equations of the form (2.2) are traditionally connected with vector fields on manifolds. In this 
particular case the right hand sides of (2.2) form the components of the vector field not on M however but on the 
tangent bundle TM 

^ ^ dx 1 dx n dv 1 dv n 



L See also |chac-dyn/9403003 and patt-sol/9404001 
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First n components of the vector field (2.3) separately can be interpreted as components of the velocity vector 

< 2 - 4 > "»'^ + - + ^ 

tangent to M . Rest part of components in (2.3) do not admit such interpretation. But if the manifold M is equipped 
with Riemannian metric g%j and with the metrical connection rfj then using the components of (2.3) one can form 
the following quantities 

(2.5) F l = & + r] k v k v j 

that do behave like the components of some vector tangent to M when we change the local map on M 

Same indices on the different levels in (2.5) and everywhere below imply summation. Vector F from (2.6) with the 
components of the form (2.5) is natural to be considered as a vector of force. The analogy of (2.1) and (2.2) then 
becomes complete. The existence of the vector F lets us speak about the angle between force and velocity. It lets 
us also break the force F into two parts directed along the velocity and perpendicular to the velocity. When F = 
equations (2.2) become the equations of geodesic line. 

3. Metric on the tangent bundle and the expansion of the algebra of tensor fields. 

Let us consider again the manifold M with the Riemannian metric gtj. Vector fields on the tangent bundle TM 
have the form 

First n components of V in (3.1) like in (2.4) are interpreted as the components of the vector X tangent to M. The 
whole set of components of (3.1) is transformed as follows 

'«> *-jsr* w ' = ^ xl + a £ w ' 

when one change the local variables x 1 , . . . ,x n for x 1 , . . . ,x n . The components of metric connection T k j under the 
same action are transformed as follows 



dx p dx" d. 



X" 



^ r P q dxidx k dx sFj?k 0x3dx k 

Comparing (3.2) and (3.3) we conclude that the following quantities Z l produced from the components of (3.1) 

(3.4) z i = w i + r l jk v k x j 

are transformed like the components of tangent vector to M under the change of local map on M. For the vector (2.3) 
the corresponding vector (3.4) coincides with (2.5). So the vector V tangent to TM can be replaced by two vectors 
X and Z tangent to M. This gives rise to the pair of linear maps ir and p where 7r is canonical projection from the 
bundle TM to the base manifold M 

(3.5) X = tt(V) Z = p(V) 

The relationship (3.4) lets us introduce the Riemannian metric on the tangent bundle TM by forming the following 
quadratic form on the set of vectors (3.1) 

(3.6) g(V, V) = g(7r(V), n(V j) + g(p(V), p(V)) 
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In terms of differentials of local coordinates the metric (3.6) is written as follows 

(3.7) g = { 9ij + rf r v r g pq r] s v s )dx l dxi + 2{g iq rJ s v s )dx i dv j + g ij dv i dv j 

Metric tensor for (3.7) is determined by the initial metric tensor g^ of M. Its matrix has the natural block structure 



(3.8) ~ 9ij 




Upper left block corresponds to the coordinates on the base while the lower right block of (3.8) corresponds to the 
coordinates on the stalk. Explicit form of metric tensor makes possible the computation of the components of metric 
connection for it but we do not need them in what follows. 

Much more interesting objects are the images of vector fields on TM under the action of maps 7r and p from (3.5). 
In general they could not be treated as vector fields on M . They are the vector-valued functions whose argument 
is the point of TM and whose value is the vector tangent to M at the point being the projection of the argument. 
Such function composes the algebra over the ring of scalar functions on TM. This algebra can easily be expanded 
up to the algebra of tensor-valued functions on M with the argument from TM. It is natural to call it the expanded 
algebra of tensor fields on M. The first example of the element of such algebra is the force field (2.6) for Newtonian 
dynamical system. This field has the real physical meaning when the manifold M is the configuration space for some 
real mechanical system restricted by inner bounds. 

Expanded algebra of tensor fields on M is equipped with the natural operations of tensor product and contraction. 
Presence of metric on M adds two operations of covariant differentiation. First of them is the differentiation by 
velocity or the velocity gradient. For scalar, vectorial and covector fields it is defined by the following expressions 

(3.9) Vnp = — ViX = — — V;A„ 



dv l dv l dv l 

Second is the covariant differentiation by coordinate or the space gradient. It is the modification of the ordinary 
covariant differentiation. The result of its application to the scalar, vectorial and covector fields is as follows 

xV dx* lk dvP 

(3.10) va m = -g- + i%x> rf k a —v* 

_ 0X m dX m ■ 

VlXm F ^ Xp ~ F ^~d^ V 

For other tensor fields the action of (3.9) and (3.10) is continued according to the condition of concordance with 
the operations of tensor product and contraction. Ordinary tensor fields are the part of expanded algebra. Velocity 
gradient for them is always zero while the space gradient coincides with ordinary covariant derivative. Particularly 

(3.11) Vfefify = V k gij = 

The relationships (3.11) express the compatibility of metric and metrical connection in terms of the above covariant 
derivatives. 

One more element of the expanded algebra is the vector field of the velocity (2.4). For the velocity and space 
gradients of it we have 

(3.12) Vfcw 4 - V k v l = 4 
Scalar field of modulus of velocity is defined by v according to the following formula 



(3.13) 



v 2 = |v| 2 = g ij v i v j 
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For the gradients of the scalar field defined by (3.13) we have 

(3.14) V k v = V k v = N k = g kq Ni 

The quantities N q in (3.14) are the components of the unit vector field N directed along the velocity v = uN. 
Gradients for this vector field are the following 

(3.15) V fc N' = VkN^v-HSl-NkN 1 ) 

They are calculated on the base of (3.12) and (3.14). Components of the matrix P£ = 5\ — N k N l in (3.15) are the 
components of operator valued field P of normal projectors on a hyperplane that is perpendicular to v. Covariant 
derivatives for P itself are 

(3.16) V fc Pj = V fc Pj = -v-\g lq P«N* + N,Pl) 

Along with P we define the additional projector- valued field with the components Qj = NjN 1 . For it we have 
P + Q = 1. Covariant derivatives of Q are easily calculated from (3.16) 

(3.17) V fc Qj - V fe Qj = v-^g^N 1 + N^) 

In addition to the properties (3.16) and (3.17) we can see that projector-valued fields P and Q are symmetric respective 
to the metric gij on M . 

Let the functions x 1 (t), . . . , x n {t) define the parametric curve on M in local map. The derivatives d t x % define the 
tangent vector to M (it may be treated as velocity vector of a point moving along this curve). Suppose that at each 
point of this curve we have the tangent vector U to M (but possibly not tangent to the curve) with the components 
U l {t). In other words we have the vector- valued function on the curve. Let us produce another vector- valued function 
on that curve according to the formula 

(3.18) v t u l = dtir + r] k u 3 d t x k 

Formula (3.18) defines the covariant derivative of the vector-valued function on the curve by the parameter t of it. It 
is well known that such derivative is zero by parallel displacement of the vector along the curve. 

Let's add the functions w 1 (t), • • ■ , v n (t) to x 1 (i), . . . , x n (t) and consider all them as the curve lifted from M to TM. 
Such lifting is called natural if v l (t) = dtx l {t). However here we consider arbitrary (may be not natural) liftings. Let us 
define some vector field U of expanded algebra in some neighborhood of the lifted curve. Substituting x 1 ^), . . . , x n (t), 
w 1 (t), . . . ,v n (t) for its argument we obtain the vector-valued function on the former curve. For this function we have 

(3.19) V t W = V k Wd t x k + V k U l V t v k 

Formula (3.19) approves the names velocity and space gradients for the covariant derivatives (3.9) and (3.10). It is 
easily modified for the case of scalar, covectorial and all other types of tensor fields of expanded algebra. 

Let us find the relation between the ordinary covariant derivatives on TM and the modified covariant derivatives 
(3.9) and (3.10) on the manifold M itself. In order to do it consider the pair of vector fields X and Y on the tangent 
bundle TM. For the projections ir and p from (3.5) applied to the commutator of these vector fields we derive 

4[X, Y]) = V t(X )7t(Y) - V T(Y) ^(X)+ 
1 ■' +V p( x)7r(Y)-V w(Y) ^(X) 

P([X, Y]) = V w( x)P(Y) - V 7r(Y) p(X)+ 
{A ' 2l> + V p( x)P(Y) - V p(Y)P (X) - R(n(X),n(Y))v 



Here R(A, B) is the operator-valued skew-symmetric bilinear form defined by the curvature tensor of M. Formulae 
(3.20) and (3.21) are proved by the direct calculations in the coordinates. 
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Let Z be one more vector field on TM and ip be the scalar field on TM which also may be treated as a scalar field 
of the expanded algebra on M. Then 

(3.22) V x ^ = a x ^ = V T(X )^ + V p(X )^ 

Left hand side of (3.22) contain the ordinary covariant derivative on TM which coincides for the scalar function with 
the derivative along the vector X. Covariant derivatives in the right hand sides of (3.20), (3.21) and (3.22) are that 
of (3.9) and (3.10). For to calculate the covariant derivatives on TM we use the following formula from [6] 

2g(V x Y, Z) = d x .g(Y, Z) + c^X, Z) - 5 z .g(X, Y)+ 
(3 - 23) + g([Z, X], Y) + ~g([Z, Y], X) + ~g([X, Y], Z) 

Covariant derivative on TM in (3.23) is ordinary one. Now let us use the relationship 

(3.24) g(X,Y) = 5 (tt(X), tt(Y)) + fl (p(X),p(Y)) 

The relationship (3.24) is derived directly from (3.6). As the result of substitution of (3.24) and (3.20), (3.21), (3.22) 
into the identity (3.23) we get two formulae 

tt(VxY) = V 7r(X )7r(Y) + V„ ( x)7r(Y)- 

- l(fl(p(Y,v)7r(X) +i2(p(X),v)7r(Y)) 

p(VxY) = V, (x) p(Y) + V„ ( x)P(Y)- 
(3.26) i 

--i2(7r(X),7r(Y))v 

Covariant derivatives in left hand sides of (3.25) and (3.26) are ordinary ones like in (3.23). Covariant derivatives in 
right hand sides of (3.25) and (3.26) are expanded ones. They should be treated as in (3.9) and (3.10). 

4. Variations of trajectories and the equations of weak normality. 

Let us consider the Newtonian dynamical system (2.2). In the second of the equations (2.2) we change the ordinary 
derivative of the velocity by its covariant derivative according to the formula (3.18). This gives us the equations 

(4.1) d t x l = v l V t «* = F i 

containing the force vector of (2.5). The equations (4.1) are more natural form for the Newton's second law on the 
manifold M. 

For the Newtonian dynamical system (4.1) now we consider the Cauchy problem with the following initial data 

(4.2) x\ =0 = x\s) ^ 1 1=0 

depending on some parameter s. Because of (4.2) the trajectories of the dynamical system also depend on s. In local 
coordinates they are given by the functions x % = x % {t, s). Their derivatives by s are the coordinates of some vector r 
tangent to M. It is the vector of variations of trajectories r l (t, s) = d s x l (t, s). The time derivative of r is the vector 
of variations of velocities 

(4.3) v t r* = d t T l + r; k r^v k 

We shall use (4.3) to obtain the equations for the vector r from the equations (4.1) of the dynamical system itself. In 
order to do it let us differentiate (4.1) by s and combine the result with (4.3) and the time derivative of (4.3). Then 
we obtain 

(4.4) V t V t T l + B i jkq v i v k T q = V 9 FV« + V 9 F 4 V t r« 
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Here R) qk is the curvature tensor for M. So the variation vector r satisfies the ordinary differential equations of the 
second order (4.4). 

The normal shift condition according to [1] and [2] consists in the orthogonality of r and the vector of velocity v. 
It is convenient to express v via the unit vector N from (3.14) and (3.15). According to the method developed in [1] 
and [2] we introduce the function ip as a scalar product 

(4.5) ip =( T ,N)=T i N i =g ij T i Ni 
For its time derivative by differentiating (4.5) we obtain 

d t v = g t] V t T l Ni +g t]T l V t Ni 

Because of (3.11) and (3.19) the metric tensor makes no contribution by covariant differentiation. To calculate V*iV J ' 
we use the relationships (3.15) and (3.19). As a result we have 

d t <fi = <feV t rW + v- 1 g l3 r l PlF k 

Because P is the symmetric operator field we may rewrite this expression in the following form 

(4.6) dt<p = g t3 V t T l N3 + v- 1 g ij F i piT« 
Differentiating (4.6) by t we obtain the second derivative for ip 

(4.7) d tt V = dtisijVt^Ni) + dtiv^g^FPy) 
Taking into account (3.15) and (3.19) we may write the first summand in (4.7) as 

(4.8) d t (g l0 V t T l N3) = War'Ni + v^F^V^ 

The first summand in (4.8) in turn is transformed by use of the equation (4.4) for the vector of variation of trajectory 

(4.9) Vut'N, = Ni\7 q F i T q + NtV.F'VtT* - N i E> agP T*v a v fl 

Let us insert the projectors P and Q into all terms in the right hand side of (4.9). In order to do it we use the 
decomposition of identical operator as 1 = P + Q. For the first summand we have 

(4.10) N i V q F i T q = \7 l F k N k P^r q + \7 q F k N q N kT I N l 
For the second summand in (4.9) by the same way we derive 

(4.11) NtVgF'VtT* = V,F k N k P l q V t T q + V q F k N q N k V t T l N t 

The third summand in (4.9) vanishes because the vectors N and v are collinear v = |v|N. Indeed 

(4.12) N i E* agfl T*v a v fl = \w\ 2 R iaqP N i N a N^T q 

Curvature tensor is skew-symmetric in i and a. Therefore the result of contraction in (4.12) is zero. 
Now let us transform the second summand in (4.7) using the relationships (3.14) and (3.19) 

(4.13) d(v- 1 g ij F i Pi T k ) = -v- 2 F k N k F l P l q r q + v^F^V^ + v- 1 W t {F i P^T q 

Then let us transform the last summand in (4.13) with the help of (3.16) and (3.19) 

v- 1 Wt(F i P^)T q = V k F l N k P'T q + v-^ k F i F k P i q T q - 
4.14 q q . q 

- v- 2 F k N k F t P t T q - v- 2 F 1 P 3 k F k T I N l 
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As a result the second summand in (4.7) may be obtained by substituting (4.14) into (4.13) and the first one may be 
obtained by substitution of (4.10) and (4.11) into (4.9) followed by substitution of (4.9) into (4.8). By analyzing the 
obtained formulae (4.5), (4.6) and (4.7) for ip, d t ip and d t t<p we find that all they are the linear functionals with respect 
to vectors boldsymbolr and V^r forming the phase space for the integral trajectories of the system of differential 
equations (4.4). Using <p, dtp and duf we construct another functional L of the form 

(4.15) du<p ~ PdtV -Q<P = L 
The coefficients for L in the formula (4.15) we define as 

P = V q F k N q N k Q = V q F k N q N k - V - 2 P q F q F k 

Then the functional L itself may be written in the form 

(4.16) L = ZiPfrtT" + Qpy 

The coefficients & and Q in (4.16) arc defined by the above calculations. They are the following 

(4.17) & = V l F k N k + 2v~ 1 F i 

, N Ci = (V^ + V k F, - 2v- 2 F l F k )N k + 

4-18 . - . - 

+ v-\y k FiF k - V k FiN k N q Fi) 

For the dynamical system (4.1) to accept the normal shift on M (see [1] and [2]) the functional L should identically 
vanish. This condition in [1] and [2] is called the condition of weak normality. In the present situation it gives us the 
following equations for & and d from (4.17) and (4.18) 

(4.20) til* = QP q = 
The equations (4.20) may be rewritten in the following explicit form 

(4.21) (v^Fi + V t {F k N k ))P; = 

(V t F k + V fc F, - 2v- 2 F t F k )N k Pl+ 

l (V k F % F k - V k F r N k N r F i )Pi = 

The total number of the equations (4.21) and (4.22) is 2n. It coincides with the twiced dimension of the manifold M. 
However because of presence of projector matrices P q in them they are not independent. The number of independent 
equations in the system of (4.21) and (4.22) is 2n — 1 which is in concordance with the results of [1] and [2]. 

The equations (4.21) and (4.22) are the covariant form of the equations of weak normality on an arbitrary Rie- 
mannian manifold. In Euclidean case M = R" they were derived in [1] and [2] by use of spherical coordinates in 
the space of velocities. The question of introducing the proper spherical coordinates here in the general situation is 
interesting but it is the subject for separate paper. Analyzing the equations (4.21) one can see that if the vector of 
force is decomposed into two parts first being along the velocity and second being perpendicular to it then the second 
part is defined by the first one. This fact was observed in [1] and [2]. It remains true for the general non-Euclidean 
situation. 

This paper is written under the financial support of two of authors (Boldin A.Yu. and Sharipov R.A.) by Russian 
Fund for Fundamental Researches. 
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